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ABSTRACT 

We use high spatial resolution (^7pc) CARMA observations to derive detailed properties for 8 
giant molecular clouds (GMCs) at a galactocentric radius corresponding to approximately two CO 
scale lengths, or ~ 0.5 optical radii (^5), in the Local Group spiral galaxy M33. At this radius, 
molecular gas fraction, dust-to-gas ratio and metallicity are much lower than in the inner part of M33 
or in a typical spiral galaxy. This allows us to probe the impact of environment on GMC properties 
by comparing our measurements to previous data from the inner disk of M33, the Milky Way and 
other nearby galaxies. The outer disk clouds roughly fall on the size-linewidth relation defined by 
extragalactic GMCs, but are slightly displaced from the luminosity-virial mass relation in the sense of 
having high CO luminosity compared to the inferred virial mass. This implies a different CO-to-B^ 
conversion factor, which is on average a factor of two lower than the inner disk and the extragalactic 
average. We attribute this to significantly higher measured brightness temperatures of the outer disk 
clouds compared to the ancillary sample of GMCs, which is likely an effect of enhanced radiation levels 
due to massive star formation in the vicinity of our target field. Apart from brightness temperature, 
the properties we determine for the outer disk GMCs in M33 do not differ significantly from those 
of our comparison sample. In particular, the combined sample of inner and outer disk M33 clouds 
covers roughly the same range in size, linewidth, virial mass and CO luminosity than the sample 
of Milky Way GMCs. When compared to the inner disk clouds in M33, however, we find even the 
brightest outer disk clouds to be smaller than most of their inner disk counterparts. This may be due 
to incomplete sampling or a potentially steeper cloud mass function at larger radii. 

Subject headings: galaxies: individual (M33) — galaxies: ISM — ISM: clouds 



1. INTRODUCTION 

Star formation in normal, non-starburst galaxies is ob- 
served to occur predominantly in giant molecular clouds 
(GMCs), which are gravitationally bound structures of 
molecular gas. In order to understand star formation 
in galaxies, it is thus important to gain a thorough un- 
derstanding of GMC properties and how they depend 
on environment within galaxies. In this paper we study 
the properties of GMCs in the outer disk of the local 
spiral galaxy M33 and study the impact of environment 
on GMC properties by comparing our measurements to 
those from the inner disk of M33 and from other local 
galaxies. 

When GMCs were systematically observed in the 
Milky Way, their observable properties were f ound to 
obey a number of scaling relations ('Larson laws '. iLar son! 
1981): Size, line width, and CO luminosity are correlated 
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in Ga lactic clouds (e.g.. lSolomon et aTlll987l : rHever et al.l 
2009). The original form of the size-line width relation 
(R oc a ' 5 ) in combination with the observation that 
molecular clouds are approximately in virial equilibrium 
also implies that their molecular gas s urface densities 
are ro ughly constant. In the Milky Way, Solomon et al. 
(1987) find E H2 ~ 17OM pc" 2 from ' virial mass esti- 
mates using 12 CO emission, though recently iHever et al.l 
(2009) found a much lower value of Sh2 ~ 50M Q pc -2 
with a scatter of 23M Q pc~ 2 using more optically thin 
13 CO emission and assuming local thermodynamic equi- 
librium. The authors of the latter study point out, how- 
ever, that this value could be as high as ~ 120M© pc~ 2 
due to abundance variations in outer cloud envelopes. 

With the advent of millimeter- wave interferometry, de- 
tailed studies of GMC properties became feasible for 
Local Group and relatively nearby galaxies. This in - 
cludes observations in, e.g., M31 flVogel et al.l [19871) . 
M33 (IWilson fc Scovilld WM lEngargiola et al.l I2003h 
and other nearby ga l axies (see IWalter et al.ll2001l 120021 : 
IBolatto et all 120081: IBlitz et al.l 120071 and references 
therein). These studies did not find large differences re- 
garding GMC properties compared to the Milky Way. 

The currently probed range in environmental parame- 
ter space is quite limited, however. Variations of cloud 
properties with environment might be expected though, 
as these properties might depend on, e.g., metallic- 
ity, interstellar r adiation field or dust abundance (e.g. 
lElmegreenl I1989D . Few studies have addressed this is- 
sue in more extreme, in particular lower metallicity, en- 
vironments and found only minor deviations (most no- 
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tably in the Small Magellanic Cloud, SMC) from the pre- 
vious l y established scaling rel a tions (e.g.. IWalter et al.l 
[20011 |200l iLerov et al.l l200l iBolatto et al.l 120081 re- 
ferred to as B08 in the following). 

After extending parameter space by going to more 
metal poor galaxies, one of the next logical steps is 
to look for varying cloud properties within the same 
galaxy. Deep single dish observations revealed CO emis- 
sion far out in or even beyond the optical disks o f 
galaxies fe.g. JBraine J^Herpin 2001 lBraine"et~aT][200l . 
including the Local Group spiral M33 ([Gardan et al.l 
l2007t iBrahTe et al.ll201Cft IGratier et al.ll2010fl . These ob- 
servations probe molecular gas properties in a regime 
where conditions in the interstellar medium and thus 
CMC properties are expected to change, because of, e.g., 
changing metallicities or dust abundance . The proximity 
of M33 (D«840 kpc, iGalleti et ail 120041 ) and the avail- 
ability of comprehensive, homogeneous measurements 
and consistently determin ed GMC properties across the 
central star forming disk (jRosolowskv et al.l 120031 120071) 
makes this galaxy an ideal target for high resolution CO 
observations to probe GMC properties in the outer disk. 

In this paper we present CO observations from the 
Combined Array for Research in Millimeter Wave As- 
tronomy (CARMA) of 8 GMCs. The galactocentric dis- 
tance of these clouds corresponds to approximately two 
scale lengths of the molecular gas, i.e., the length over 
which the azimuthally averaged CO emission declines by 
a factor e 2 , or ~ 0.5r25, where 125 is the isophotal radius 
corresponding to 25 B-band magnitudes per square arc- 
second. Figure [Tj illustrates the location of our targeted 
field and shows the CO peak intensity map we obtain. 
We introduce our observations which we use to measure 
cloud properties for 8 GMCs in Section [2] and discuss the 
expected environmental variations (molecular gas frac- 
tion, dust-to-gas ratio, metallicity) within M33 with re- 
spect to our target field in Section [3] We compare our 
measurements to those from the inner disk of M33, other 
nearby galaxies and the Milky Way in Section [4] 

2. DATA & METHODOLOGY 

Our CARMA observations were targeted towards one 
of the outermost detected mol ecular gas complexes from 
the single dish observations of iGardan et al.l (120071 com- 
pare Figure [1]). We note that recently iBraine et all 
(2010) have obtained deeper single dish data and de- 
tected CO emission even further out in the disk of M33. 
We observed the target position in the CO(1-0) transi- 
tion using a 19-point mosaic in D configuration between 
July and August 2008 and in C configuration between 
October and November 2009 under mostly good 3mm 
weather conditions. The total observing time was ~37hr 
and ^30 hr in D and C configuration, respectively. The 
synthesized beam size is ~ 1.7", corresponding to ~ 7pc 
at the distance of M33. We tuned the 3mm receivers to 
the doppler-shifted CO(1-0) rest frequency in the upper 
sideband and placed two ~ 31 Mhz wide bands across the 
line. This yields an effective bandwidth of ^155 km s -1 
and a velocity resolution <5v c hannci = 1.27kms . For 
phase and amplitude calibration, 0205+322 was ob- 
served every 20 minutes. Fluxes were bootstrapped from 
Uranus, Neptune or MWC349 and 3C454.3 or 3C84 were 
observed for bandpass calibration. Radio and optical 
pointing was performed every 2 hr. We estimate the re- 



sulting calibration to be accurate within ~ 15%. Data 
reduction was performed in MIRIAD. Figure Q] shows the 
target field with respect to the young stars (left panel) 
and the CO(1-0) peak intensity map from our observa- 
tions (right panel). 

We identify GMCs a nd measure their properties using 
the CPROPS package (|Rosolowskv fc Lerovll2006l) . We 
define GMCs as contiguous regions of high signal-to- noise 
emission. No further decomposition appears necessary in 
this datasct. CPROPS uses moment methods to measure 
radius, line width, luminosity, and peak temperature for 
each cloud. It estimates the associated uncertainties us- 
ing boostrapping techniques and attempts to correct each 
measurement for biases due to finite sensitivity and lim- 
ited spectral and angular resolution. A key point for this 
study is that B08 have used exactly the same approach to 
measure GMC properties from a large sample of nearby 
galaxies, including the inner disk of M33. This gives us 
a well-controlled point of comparison. 

Table [1] lists the derived properties for the 8 GMCs in 
our field. We note that more faint CO features are visible 
in Figure [TJ but the signal-to- noise ratio in these features 
is too low for reliable cloud property estimates. 

3. ENVIRONMENT 

Previous interferometric observations of individual 
GMCs in M33 focussed on GMC properties i n the 
inner dislfl of M33 (e.g ., lEngargiola et~aTl 12001 
iRosolowskv et al.ll2003L 120071 ) . In the following we assess 
the distinguishing characteristics of the ISM at larger 
radii in M33 and how this may impact GMC properties. 

Numerous studies have compared the radial behavior 
of atomic gas, molecul ar gas and various star forma- 
tion tracers in M33 (e.g . . lEngargiola et al.ll2003t iCorbellil 
12001: IHever et al.ll200l IGardan et al.ll2007D . While the 
H I is found to remain relatively constant across the disk 
of M33, these studies derive a range of scale lengths be- 
tween ~ 1.4 — 2.5 kpc for the CO emission. We note 
that we adopt a value of 2 kpc for the CO scale length 
throughout this paper. From the different behavior of 
H I and CO emission, it is clear that the molecular 
gas fraction E h 2 /£hi decrea s es as a function of radius. 
Most recently, IGratier et al.l (|2010[) measured this frac- 
tion and find that while the gas is mostly molecular in 
the center, the molecular gas fraction decreases to only 
i~ 10 — 20% at our target position at r« 4 kpc (adjusted 
to our adopted Galactic CO-to-H2 conversion factor of 
X co = 2 x 10 20 cm- 2 (Kkms- 1 )- 1 ). 

The stellar surface density, £ s tar, in M33 is low com - 
pared to other, larger spirals (e.g., IBraine "eTaI]|20T0l) . 
As S s tar is important for the hydrostatic gas pressure, 
which was found to correlate wel l with the molecula r 
gas fraction in nearby galaxies (e.g-. IWong fc B litz 2002), 
this may be an explanation for the unusually low H2 frac- 
tion. Another factor that may impact the H2 fraction is 
the amount of dust in the ISM, as dust grains serve as 
sites of H2 formation and s hield the molecular hydro- 
gen from photodissociation. IGratier et al.l (|2010t ) derive 
24, 70 and 160 /im scale lengths of ~ 1.4, ~ 1.5 and 
~ 1.8 kpc, respectively, all of which are quite similar to 

8 We compare our measurements to the M33 GMC sample from 
IBolatto et al.l 120081) . which has a median galactocentric distance 
of less than one CO scalelength. 
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Fig. 1. — Left: The location of the GMC complex analyzed in this paper (black rectangle; the size matches the peak intensity map in 
the right panel) relative to the distribution of young stars (DSS B-band image) in M33. The galactocentric distance of our target field 
corresponds to approximately two CO scale lengths (the white ellipse indicates one CO scale length or r =2kpc). It contai ns one of the 
furthe st, large GMC complexes in the disk of M33 that is detected with high signal-to-noise in deep single dish observations by Card an et aO 
120071) . Right: 12 CO(1-0) peak intensity map from CARMA of our target field. The size of the synthesized beam of ~ 1.7" is shown in 
the lower left corner. CPROPS identifies 8 GMCs in this field (highlighted by contours), for which, at ~ 7pc resolution and with peak 
signal-to-noise ratios > 10, we can derive accurate cloud properties. The number next to each cloud identifies them in Table [T] where we 
present the derived properties. 



TABLE 1 
GMC Properties 



Number a 


R.A. 


Dec. 


VJsr 


R b 




M ium b 


M vir b 


T B C 


SGMC d 




[J2000] 


[J2000] 


[kms -1 ] 


[PC] 


[kms -1 ] 


[1O 3 M ] 


[1O 3 M ] 


[K] 


[M pc" 2 ] 


1 


01 33 59.0 


30 55 32.5 


-254.6 


10.1±1.9 


1.6±0.5 


47±4 


27±18 


10.4±0.6 


84±64 


2 


01 34 01.9 


30 54 34.2 


-257.4 


18.3±1.4 


1.6±0.3 


105±7 


46±20 


8.7±0.6 


44±20 


3 


01 34 02.7 


30 53 48.3 


-258.1 


12.3±1.4 


1.5±0.3 


76±9 


30±14 


11.0±0.6 


64±32 


4 


01 34 03.0 


30 55 05.9 


-256.5 


14.8±2.1 


2.1±0.5 


62±9 


70±39 


6.5±0.6 


102±64 


5 


01 34 04.2 


30 55 06.5 


-256.8 


22.9±1.6 


2.0±0.3 


233±9 


91±29 


9.9±0.6 


55±19 


6 


01 34 04.2 


30 55 19.3 


-256.6 


13.4±1.8 


2.7±0.5 


60±8 


104±48 


5.7±0.6 


185±100 


7 


01 34 04.7 


30 54 26.4 


-255.2 


16.1±1.3 


2.3±0.4 


91±9 


90±36 


6.1±0.6 


111±47 


8 


01 34 04.9 


30 54 39.2 


-265.2 


17.0±1.5 


0.9±0.1 


84±5 


14±2 


7.0±0.6 


15±4 



a The numbers refer to the individual clouds in the right panel of Figure [T] 
b Uncertainties are estimated from bootstrapping in CPROPS. 

c The uncertainty on Tb is estimated from the RMS scatter of the noise in the data cube. 

d Sgmc i s derived from the virial masses and radii in this table. Uncertainties are estimated from error propagation. 



the CO scale length in M33. In combination with the flat 
H I profile, this indicates 1) that the dust-to-gas ratio is 
significantly lower (~ 25%) at r«4kpc compared to the 
inner disk and 2) that the dust-to-gas ratio may play an 
important role in setting the molecular gas fraction. 

Because the dust-to-gas ratio also scales with metallic- 
ity, one might also expect a lower metallicity in our tar- 
get region. Even though earlier measurements suggested 
a quite significant m etallicity gradient across M33 (up 
to ~ -0.1 dexkpc -1 , iVilchez et all 119881: iGarnett et all 
1997), recent measurements provide evidence for a much 
shallower gradient of < — 0.03 de xkpc -1 ( Crocket t et al.l 
l2006HRosolowskv fc Sim on 2008). Compared to the cen- 
ter, this implies a metallicity that is lower by ~ 30% at 
r « 4kpc, i.e., Z» 8.24, ad o pting a central value of 8.36 
from iRosolowsky &; Simonl (|2008f ). This metallicity cor- 
responds to only ~ 30% of the Galactic average. 

A lower metallicity, i.e., lower carbon and oxygen abun- 



dance, and dust-to-gas ratio in the ISM are likely to affect 
GMC properties, because they lead to lower H 2 forma- 
tion r ates and less e ffective shielding f rom UV radiation 
(e.g.. lMcKee1ll98i lElmegreenl [19891 : iMalonev fc Blackl 
1988). Thus, one might expect GMCs to have higher 
(surface) densities in such environments and possibly 
to consist of smaller CO-bright cores embedded in CO- 
dark H 2 envelope s, leading to a higher CO-to-H 2 conver- 
sion factor (e.g iMadden et"all 119971: iLerov et all 120071: 
lElmegreenl fl989f ). Our new measurements in combina- 
tion with ancillary data allow us to test these expecta- 
tions in M33 and to extend such studies into a barely 
probed regime of ISM parameter space. 

4. GMC PROPERTIES AND COMPARISON TO OTHER 
MEASUREMENTS 

Figure [5] shows spectra for the 8 identified GMCs. 
These clouds have relatively narrow spectra and the de- 
rived line widths are notably smaller than those from 
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Fig. 2. — Spectra of integrated emission within the isophotal contours in the right panel of Figure [T] for each cloud in each channel. 
All clouds are clearly detected in several channels. The cloud numbers correspond to the numbers in Table [T] listing the derived cloud 
properties. 

clouds from the inner disk of M33 (see discussion for Fig- 
ur e [H below) . This agr ees qualitatively with the finding 
of iBraine et al.l (|2010f) , though their spatial resolution 
is much lower, which may lead to overestimating line 
widths by averaging over multiple clouds. Most of our 
clouds have peak flux densities between ^1 — 3 Jy, with 
the exception of cloud 5 (~ 5 Jy), which is the biggest 
and most massive GMC in our sample (though note the 
similarly small inferred velocity dispersion of ~ 2kms _1 
compared to the other clouds in Table [I]) . 

B08 have assembled and partly reanalyzed CO obser- 
vations of 3 disk galaxies (the Milky Way and the Lo- 
cal Group spirals M31 and M33) and 11 Local Group 
and nearby dwarf galaxies. Because in this paper we 
follow the same methodology, we can directly compare 
the properties of the B08 compilation of GMCs to our 
outer disk GMCs in M33. Comparing the peak bright- 
ness temperatures, we find a mean value for the new 
GMCs that is significantly higher than in the B08 en- 
semble: TB,ncw = 8.2 K with a ler rms scatter of 3.3 K 
and T b ,bo8 =2.3K with a scatter of 1.2K. For M33 
specifically, all of the outer disk clouds have higher Tb 
than all but one of the inner disk clouds. Tb depends 
mainly on the beam filling fraction and the excitation 
temperature (it also depends on optical depth, but the 
12 CO(1-0) emission is most likely optically thick in all 
cases). In B08, typical brightness temperatures within 
individual galaxies rarely exceed 3.5 K, even for cases 
where GMCs were observed at high spatial resolution 
(comparable to this study) so that GMCs were comfort- 
ably resolved. This also includes the previous, inner disk 
M33 observations. At least for these higher resolution 
data, the equally high beam filling fractions > 1 argue 
for higher kinetic temperatures in the molecular gas (un- 
der the assumption of uniform brightness distributions 
within GMCs) to explain the higher peak brightness tem- 
peratures we measure. A possible explanation might be 
the potentially elevated radiation field due to massive 
star formation in our target field. 
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Fig. 3. — Comparison of 24/^m (red), far UV (blue) and Ha emis- 
sion (green) in our target region. The full sensitivity field-of-view 
is indicated with a white line and integrated intensity contours of 
our CO observations are shown in cyan. The contours are running 
from 0.25-1.2 Jy beam - km s in steps of 0.2 Jy beam" 1 kms . 
The size of the synthesized beam is shown in the lower left cor- 
ner. The cloud numbers are identical to those in Figure [T] and 
refer to the cloud properties in Table [l] The figure shows several 
regions of intense 24/im and Ha emission (both indicating current 
star formation) in the immediate vicinity of the detected GMCs. 

We illustrate this in Figure where we show 

24/iin (red), far U V (blue) and Ha emission 

(jHoopes fc Walterbosl 120001 green) in our target 
region. Overlaid are the CO integrated intensity 
contours in cyan and our full sensitivity field-of-view 
(white line). From this figure it becomes clear that 
our GMCs are in the immediate vicinity of several 
star forming regions. In particular the two clouds with 
the highest peak brightness temperatures, 1 and 3, 
are directly associated with Ha and 24fim sources. 
The combination of presumably low dust-to-gas ratios, 
low gas columns and intense star formation supports 
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Fig. 4. — Size-line width (left panel) and CO luminosity-virial mass (right panel) plots f or the new M33 clouds ( black points, this study) 
as well as for a number of other datasets for com parison: inner disk M33 (gray points, Rosolowsky ct al. 2003, BOS) and SMC clouds 
(triangles, B08), Milky Way GMCs (black crosses. ISoIomon et al. 1987; Hcvcr ct al. 2009, B08) and a fit to an ensemble of extragalactic 
CMCs (dashed line, B08). The new clouds extend the distribution of the previously known clouds in the inner disk of M33 toward smaller 
sizes and line widths as well as toward smaller virial masses and lower CO luminosities. The combined M33 inner and outer disk sample 
covers about the same phase space in both plots as the Milky Way clouds. 



our interpretation that the molecular gas is heated to 
high temperatures leading to the unusually high peak 
brightness temperatures we measure. 

In FigureS]we analyze scaling relations for our GMCs: 
we plot the size-linewidth relation in the left panel and 
the relation between luminosity and virial mass in the 
right panel. In both panels, the dashed line shows the 
best fit to the cloud ensemble from B08 and the gray 
points represent GMCs from the inner disk of M33 from 
I Rosolowsky et al.l (|2003l ). but as reanalyzed in B08 using 
CPROPS. We also add clouds from the Small Magellanic 
Cloud (SMC; triangles), the most metal-poor system in 
the B08 sample (Zsmc ~ 0.2 Z mw)- Black crosse s rep- 
resent Milky Way GMCs fromjHever et al.l (12009T) (who 
reanalyzed the ISoIomon et al.l ( 19871) sample) , which in- 
cludes a number of small clouds quite similar in size to 
our sample. We reiterate that they use a somewhat dif- 
ferent methodology and a different tracer ( 13 CO emis- 
sion) compared to the other studies (see Section [TJ. The 
black points show the new M33 outer disk GMCs from 
Table [TJ The error bars show estimates of the uncertain- 
ties from bootstrapping in CPROPS. 

The left panel shows that 1) the new GMCs roughly fall 
on the B08 fit, but that 2) the data points are offset com- 
pared to the majority of inner disk clouds in M33. Due 
to limited sensitivity and resolution in the extragalac- 
tic measurements, however, small GMCs are not always 
readily observable. For the Milky Way, however, where 
observations do have sufficient resolution and sensitivity, 
the GMCs overlap the combined inner and outer disk 
M33 points, indicating that the outer disk GMCs do not 
show dramatically different properties compared to what 
has been measured in the Galaxy. Nonetheless, most of 
the clouds we observe at ~ 0.5r 2 5 are smaller and tend to 
show lower velocity dispersions than most of the clouds 
in the inner disk of M33. This might hint at either a 
steeper cloud mass function or incomplete sampling of 



the latter. 

We compute surface densities for our sample via 
M v i r (7rr 2 ) -1 and derive a mean surface density of 
82 M pc" 2 with a la rms scatter of 43 M pc" 2 . This 
value is well within the range of values quoted in the lit- 
erature, 5OM pc" 2 < S H 2 < 170M Q pc" 2 , including 
Galactic clouds (compare Section [JJ or the extragalactic 
sample in B08. Thus, we find quite typical surface den- 
sities rather than much larger values as might have been 
expected (Section [3]). 

The CO luminosity-virial mass plot in the right panel 
shows that the outer disk M33 GMCs seem to be an ex- 
tension of the distribution of inner disk clouds toward 
lower CO luminosities and virial masses. As in the left 
panel, the combined M33 sample overlaps the Milky Way 
ensemble, indicating a similar range of virial and lumi- 
nous GMC masses in both spiral galaxies. The new M33 
clouds appear to fall slightly below the ensemble fit of 
B08, however. Because the normalization (intercept) of 
the power-law fit relating virial mass to CO luminosity 
yields the CO-to-H 2 conversion factor, aco, the offset 
of our data points from the B08 fit implies a different 
average aco (assuming identical power law slopes). 

In order to assess this quantitatively, we compute 
aco from the virial and luminous masses for each 
of the inner and outer disk clouds in M33 and find 
mean values of 6.8 ± 0.8 M (K km s" 1 pc 2 )" 1 and 3.3 ± 
O.8M (Kkms" 1 pc 2 )" 1 , respectively (the quoted er- 
rors represent the la uncertainty in the mean). B08 
derive 7.6^2 6 M (K km s" 1 pc 2 )" 1 for the extragalactic 
ensemble (the errors represent the scatter in the data), 
which is in good agreement with the M33 inner disk 
value. We test the significance of the difference between 
the M33 inner and outer disk means with a Student's 
t-test, which yields a probability of ~ 1% for the null 
hypothesis that both means are equal. This implies that 
the average CO-to-H2 conversion factor for the outer disk 
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GMCs is in fact smaller (by about a factor of two) than 
the extragalactic and the M33 inner disk average. This 
finding is contrary to the expectations based on environ- 
mental conditions, which suggested a somewhat higher 
conversion factor (Section |3|) . 

One factor that could influence aco (and is particu- 
larly relevant in outer disks) is a decreasing fraction of 
CO emitting H2 in the more dust and metal-poor en- 
vironment at larger radii (see Section [3]). We therefore 
compare our measu rements to complem entary aco esti- 
mates in M33 from iLeroy et al.l (|2010l ). They estimate 
aco from dust modeling (thus tracing the entire H2 dis- 
tribution under the assumption that gas and dust are well 
mixed) and find values of ~ 6.3M Q (K km s _1 pc 2 ) -1 
and ~ 4.7M Q (K km s -1 pc 2 ) -1 for the inner and outer 
part of M33, respectively. These numbers are in good 
agreement with our virial mass estimates, which argues 
against a significant amou nt of CO-dark H2 at large radii. 
We note that because the ILeroy et"aTI ()2010l ) values are 
averages over a large area, the application of these val- 
ues to our specific region does not rule out conclusively 
unaccounted for H2. 

Taken at face value, however, the lower average con- 
version factor we measure for the outer disk clouds could 
be interpreted as a higher fractional CO abundance, i.e., 
more "CO per H2". This is highly unlikely, however, 
given the much lower dust-to-gas ratios and metallic- 
ities at larger radii in M33 (compare Section [3]). On 
the other hand, aco scales inversely with the bright- 
ness temperature : aco <x Tb _1 (|Dickman et al.l H986t 
iMalonev fc Bla^[l98a ). Thus, the more likely explana- 
tion seems to be the systematically higher Tb for the 
outer disk GMCs, which we interpreted above as higher 
kinetic gas temperatures due to elevated radiation levels 
from nearby star formation. The higher temperatures 
could lead to hig her excitation of the CO molecules (e.g., 
iWeifi et aT1l200lD . thus lowering the measured aco- 

With the resolution and sensitivity of CARMA we are 
able to measure the properties of 8 GMCs in the heav- 
ily H I-dominated outer part of M33. Despite an envi- 
ronment very distinct from a normal spiral galaxy (low 
molecular gas fraction, stellar surface density and dust- 
to-gas ratio), the clouds we observe show generally sim- 
ilar properties compared to GMCs in the Milky Way or 
other nearby galaxies. The main difference is that the 
gas appears to be hotter, with excitation temperatures 
between ~ 6 — 11 K, which is likely to be the responsi- 
ble mechanism for a lower inferred CO-to-H2 conversion 
factor. This difference appears mostly attributable to 
heating by massive star formation coincident or adjacent 
to the GMCs. 
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